Exact time correlation functions for iV classical Heisenberg spins in the 'squashed' 
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We present exact integral representations of the time-dependent spin-spin correlation functions for 
the classical Heisenberg iV-spin 'squashed' equivalent neighbor model, in which one spin is coupled 
via the Heisenberg exchange interaction with strength Ji to the other N — 1 spins, each of which 
is coupled via the Heisenberg exchange interaction with strength Ji to the remaining N — 2 spins. 
At low temperature T we find that the N spins oscillate in four modes, one of which is a central 
peak for a semi-infinite range of the values of the exchange coupling ratio. For the N = 4 case of 
four spins on a squashed tetrahedron, detailed numerical evaluations of these results are presented. 
As T — > oo, we calculate exactly the long-time asymptotic behavior of the correlation functions for 
arbitrary N, and compare our results with those obtained for three spins on an isosceles triangle. 



PACS numbers: 05.20.-y, 75.10.Hk, 75.75.+a, 05.45.-a 



I. INTRODUCTION 

Recently there has been a growing interest in the study 
of the properties of magnetic molecules. The 
defining characteristic of these substances is the presence 
of a small cluster of magnetic ions located at the center 
of each molecule and surrounded by a complicated struc- 
ture of non-magnetic chemical ligand groups. In general, 
the strength of the magnetic interaction between ions 
located in different molecules is negligible in compari- 
son to the strength of their intramolecular interactions. 
Therefore, measurements of the magnetic properties of 
macroscopic samples reflect the underlying magnetic in- 
teractions within a single molecule. 

The list of synthesized magnetic molecules has been 
constantly growing, even though most of the experimen- 
tal activity has been focused on the determination of the 
magnetic properties of a molecule containing twelve man- 
ganese ions at its core, often referred to as Mn.12. The 
theoretical tools currently used to describe the behavior 
of this relatively complicated structure are still rudimen- 
tary, and are based on a single-spin phenomenological 
Hamiltonian. 1] It is important to notice that a num- 
ber of molecular structures containing smaller numbers of 
magnetic ions have already been synthesized. For some of 
these structures it is possible to perform a more detailed 
theoretical analysis of their magnetic behavior starting 
from a many-spin Hamiltonian. 

Among the smaller clusters are a regular tetrahedron 
of Cr 3+ ions (S = |) 0*0], CV4, and a squashed tetra- 
hedron of Fe 3+ ions (S = §) 0,0, Fe 4 . For increas- 
ing values of the spin a description in terms of classical 
spins is expected to capture many of the features of the 
system 0, 0- In the present paper we provide exact 
expressions for the time-dependent spin-spin correlation 
functions for the classical Heisenberg A^-spin squashed 



equivalent neighbor model, which is the AT-spin general- 
ization of four classical Heisenberg spins on the corners of 
a squashed tetrahedron. Specific numerical results for the 
squashed tetrahedron case, N = 4, are provided. Analo- 
gous studies have recently appeared for three spins on an 
isosceles triangle and on a chain 8], for four spins on a 
square ring , and for the equivalent neighbor model of 
N classical spins |lOj . Quantum time-dependent correla- 
tion functions have been computed for a dimer |g and for 
three spins on an equilateral triangle , and for a dimer 
of classical and quantum spins in a constant magnetic 
field . The availability of time-dependent correlation 
functions is necessary to analyze neutron scattering ex- 
periments. 

In Section II we define the Hamiltonian system to be 
studied and write the corresponding partition function. 
In Section III we present the constraints upon the vari- 
ous correlation functions. In Section IV we present our 
analytic results for arbitrary N. We evaluate the long- 
time behavior of the correlation functions at infinite tem- 
perature T, and provide analytic formulae for the low-T 
modes for arbitrary N. In Section V, we present numeri- 
cal results at low T for the squashed tetrahedron, N = 4. 
Section VI contains our conclusions. A collection of in- 
termediate steps useful to the calculations is compiled in 
the Appendix. 



II. THE MODEL 

We consider N classical spins of unit magnitude, \Si\ = 
1, interacting according to the Hamiltonian 



H 



J: 



M 



M 
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where 



M = N - 1 > 2. 



(2) 



Introducing the total spin S — J^Li ^ an( i the auxil- 
iary variable Sx^m = YlfLx &i = & ~ "^JVi the Hamilto- 
nian can be written as 



1 a2 



(3) 



where we have dropped the constant energy (Jj + 
MJ2)/2, The partition function can then be calculated 
following the technique described in 0, [^. Letting 
s = |S| and x = \S\^m\i one obtains 

Z = dxV M {x) / sdsexp(-0H) (4) 
e Q f M 

= — / dxT>M(x) exp(a7X 2 ) sinh(2ax), (5) 
a Jo 



where = (fceT)- 1 , a = 0Ji/2, 7 = J2/J1, and V M {x) 
is the classical M-spin density of states, [lfj which we 
redisplayed in Eq. l|A44jl in the Appendix. 

In order to compute the time-dependent correlation 
functions, we first solve the classical equations of motion 
appropriate for the Hamiltonian, Eq. J2J), 



Sn,i-*m — JiSn,i->M x S 



(6) 



and S = 0, so that S is a constant of the motion. Fol- 
lowing the technique illustrated in || 0, , we obtain 

x [cos(si*)a: — sin(si*)y] , (7) 

where t* = Jit, s = S/s = xxy, C N = (s 2 -x 2 + l)/(2s), 
C^ M = (s 2 +x 2 - l)/(2a), ^ = 1 - C 2 ,, and A N = 

We must also consider the equations of motion for the 
Si(t), i = 1,2, ... ,M. In order to calculate the time 
correlation functions, symmetry allows us to choose just 
one of them, i = 1. We then write Sx~>m = Sx + S^mi 
and solve 

Si,2-*M = J^Sxp^tM X S + ( Ji — J2)S : 1,2— Af X 5jv- (8) 

After defining Si± = Sxx ± iSi y , we obtain, 
Si±(*) = — ^ exp(=F«st ) 



x exp {i [=ps + (1 - 7)21] t* + i0 o } 
_ 2(C 1 ^ M ±a ; ) 

x exp {i [ts - (1 - 7)x] i* - i0 o } , (9) 
■M*) = 5i z0 + A5i z0 cos [(1-7)^ + 00], (10) 



where <f>o is an arbitrary phase, and similar equations for 
the components of S2~,m- After combining these equa- 
tions with analogous ones for the components of Sx->m> 
the constants appearing in Eqs. © and (|10|l must satisfy 



SlzO — 



(AS lz0 ) 2 



Cx- 



At 



■M 



1 



1 
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(x 2 - y 2 + l) 2 
4a; 2 



(11) 
(12) 



where y = \S 2 ^m\- 

Previously, we solved these equations for the simplest 
case, M = 2, for which y = 1 is not a variable. In that 
case, the correlation functions were obtained from the 
double integrals over x and s, according to the weighting 
factors in Eq. For M > 3, however, y can vary over 
the entire range < y < M — 1. Hence, for the explicit 
evaluation of the correlation functions with M > 3, it is 
useful to rewrite the expression of the partition function 
(|4*|) in terms of a triple integral over s, x, and y, 

*M-1 /-y+1 rx+1 

Z= / T>M-i(y)dy / dx sds exp(—0H). 

Jo J\y-1\ J\x-1\ 

(13) 



III. CONSTRAINTS 

In this section we analyze the constraints upon the 
time-dependent spin-spin correlation functions 



C ij (t) = (S i (t)-S j (0)), 



(14) 



where the thermal average (...) is performed by averag- 
ing over the arbitrary phase </>o and the variables s, x, 
and y, with respect to the canonical ensemble defined by 
Eq. Ijl3|l . D ue to the symmetry of the molecule, only 
four of the N(N + l)/2 correlation functions in Eq. i|14|) 
are distinct. We write these as Cxx(t), Cx2{t), Cijv(i), and 
Cjvjv(*)- Conservation of the total spin adds a constraint, 



(s 2 ) = C NN (t)+MCxx(t) 
+M(M- l)C 12 (t). 



2MC 1N (t) 



(15) 



Finally, by writing the multispin correlation function 
(Sx-^u(t) • Sx^m{0)) in two ways, we find a second con- 
straint between two of the correlation functions, 



(sC N } = C NN {t)+MCx N (t), 



(16) 



where the constant Cm is given just below Eq. 0. The 
two remaining independent correlation functions C]\fN{t) 
and Cix(t) must then be calculated by explicitly substi- 
tuting into Eq. (|14fl the time dependences obtained in 
Section II. For Cjviv(i), this is relatively simple, as one 
can just use Eq. (JJJ) for Sp[(t), which is independent of 
y, to evaluate it. This leads to 



C NN {t) = (C N + A N cos(st*)), 



(17) 



3 



which can be evaluated using the simplified weighting 
factors present in Eq. Q. From Eq. i|16|) . this simplifi- 
cation also applies for CiAr(i). We note that Eq. l(T7|) dif- 
fers from the expression for the autocorrelation function 
in the TV-spin classical Heisenberg equivalent neighbor 
model only by the x dependence of the Hamiltonian. [TgT| 
which is irrelevant as T — ► oo. 

The challenge is to calculate Cu{t). It is useful to 
separate the expression for Cu(t) into the four integrals 
Ii(t) (i = 0,...,3), 



c 1 i(i)=X)/ i (t). 



(18) 



i=0 



The explicit triple integral representations of the /,-(£) 
valid for arbitrary T are given in the Appendix, where it 
is also shown how to reduce them to double integrals. 



IV. ANALYTIC RESULTS FOR ARBITRARY N 
A. Infinite temperature limit 

Here we present our results for the correlation func- 
tions with general N values as T — > oo. As shown in 
the Appendix, in the limit T — > oo, the triple integrals 
appearing in (|18|l can be reduced to single integrals. For 
N = 4, the relevant density of states appearing in Eq. 
(fT3")l is V%(x) — |0(a;)9(2 — x), so this reduction in the 
number of integrals is relatively simple. As T — > oo, the 
different couplings appearing in the Hamiltonian become 
irrelevant for Cjyjv(t), so that it becomes equivalent to 
that of the TV-spin equivalent- neighbor model. [Icj] 



lim Cn n (t) 



1/N + M[S N + f N (t)}, (19) 



where f N (t) ~ (t*)~ N for t* » 1. Since as T -> oo, 
(s 2 ) = N, (x 2 ) = M, and (y 2 ) = M - 1, from Eqs. JTHJ 
and l|19|l . we have 



lim C 1N (t) = 1/N - 8 N - f N (t). 



(20) 



For Cn(i) and Ci2(i), even as T — > oo, the situa- 
tion is more complicated, as the results depend crucially 
upon the values of 7 = Jij J\- As t — > 00, the time- 
dependent trigonometric functions in 1\ , I2 , and -Z3 os- 
cillate increasingly rapidly and yield vanishing contri- 
butions to lim^—tfy, Ci 1 (t), as stated by the Riemann- 
Lebesgue lemma [13J. Therefore, for arbitrary N, 



lim Clf{t)=h = {Sl z0 ), 



(21) 



We note that Iq depends upon N, and is a rather messy 
triple integral, but that as T — > 00, can be evaluated 
exactly, as shown in the Appendix. 



At infinite temperature one obtains for N 
1 

4 



where 64 
and 



lim Cait) 

T— »a© 

= -(11/180) 



35 4 



0.436345, 
(8/45) In 2 



(22) 

0.062115, HE3 



lim Clf\t) 



t— * 



0.355496, 



(23) 



the exact expression for which is given in 1A47|) in the 
Appendix. In Table I in the Appendix, we also list the 
T -> 00 values of lim^oo Cjf ^(t) for 3 < N < 11, and 
compare them with the T — > 00 values of lim^oo Cjvjv(i)- 
We note that as T — > 00, for each of these N val- 
ues, limj-foo Clf 1 ^) < limt^ooCNNit)- As T — > 00, 
limt^oo Cjvjv(t) decreases monotonically with increasing 
N to I as N — > 00. Since lim^oo C7/^ 1 (t) also de- 
creases monotonically with increasing A, and for 8 < 
N < 11, its value is less than i, it appears that this 
inequality is likely to hold for all N values. 

We now turn to the long-time asymptotic behavior of 
Cn(t) at infinite T. Following the method described in 
|jg, we first define SCij(t) = Cij(t) - lim^oo dj For 
7 = 0, the dominant behavior of limr-^oo Cn(t) is given 
by h(t), but for 7^ 7 7^ 1, it is given by 2a(i)- At long 
times, i 3> 1, where i = (1 — 7)i*, one can evaluate the 
asymptotic behavior as T — > 00 exactly. By integration 
by parts M times, we find, 



lim 5Cjf>\t) 



E(M/2 



An 



x cos[(M - 2p)t 4 



-2p) 

Mtt/2], (24) 



where 



/(y) = i + 2T 2 - 



(y 2 -i) 2 lr /j/ + i 



4y 3 



1 



(25) 



and Amp is given in the Appendix. Although the func- 
tion f(y) is non- analytic at y = 1, it can be shown that 
its derivatives do not contribute to the long-time asymp- 
totic behavior. In addition, for t* 1, one can easily 
obtain the asymptotic expression of I3 (t) , leading to 

Jim 8Cl=\t) ~ I™'" dyVM.riy) 

t*»i 



4i* Jo 

xy 3 /(y)- 

In particular, for A = 4, we obtain 



(26) 



lim 5C4 

T — foa 
t*>l 



^44 (t) 

hm SC^f ' 1 ^) 

t»l 

hm sejr°(t) 



4(r 
1 



- cos(4i*) 



(27) 



8(t) 3 



/(l)sin(t)+/(3)sra(3f) 



23 9 
30~40 ln3 



sin(i*) 
t* 



(28) 
(29) 



4 



where Eq. (j2"?|) was given previously. 0, ^| 

It is interesting to compare the present results to the 
analogous ones obtained for the isosceles triangle of spins, 
N = 3. 8] For N > 4, the infinite-T, long-time behavior 
of C\i(t) for 7 ^ 0,1 is determined by the integral hit) 
given by Eq. I|A41() . For N — 3, an additional contribu- 
tion to the infinite-T, long-time behavior of C\\(t) arises 
from I 3 (t) given by Eq. (fA"42")) . 8] For TV > 3, the cor- 
relation function Cu(i) for 7 ^ 0, 1 decays slower than 
Cn n (t) [denoted C22W in Ref. Q for N = 3], approach- 
ing its long-time asymptotic value at infinite temperature 
as (t*)~ M . In the limiting situation 7 = 0, corresponding 
for N — 3 to the three-spin chain (or 'two-pronged star') 
and for N > 4 to an M -pronged star of spins equally 
coupled to a central one, as T — > 00 and t* ^> 1, the 
correlation function is dominated by Ia(t). In this case, 
Cn(t) approaches its asymptotic limit much more slowly, 
as (i*) -1 , as shown in Eq. l(2"6"j> . 



B. Low-temperature correlation functions 

At any finite temperature, it is not possible to reduce 
the time-dependent correlation functions to a single in- 
tegral representation, even for N — 4. Since the time- 
dependence of the integrand is a simple trigonometric 
function, it is convenient to compute the Fourier trans- 
forms of the 5Cij(t), quantities which are anyhow of di- 
rect experimental relevance in neutron scattering exper- 
iments. In this case, it is then possible to express the 
Fourier transforms in terms of a single integral represen- 
tation, which then allows a precise and fast numerical 
integration. We limit our numerical work to the case of 
the squashed tetrahedron, N — 4 (M = 3). 

We define the Fourier transform as usual as 



7T 



dteyLv(iujt)8C i:j {i). (30) 



The position of the various peaks as a function of 7 
may be obtained analytically in the T — > 00 limit through 
an asymptotic evaluation of the integrals, or numerically 
by plotting the curves at large enough values of \a\ oc 
1 /T. In the Appendix, we have sketched the derivation of 
the low-temperature mode frequencies for general N, for 
both FM and AFM cases. For ferromagnetic couplings, 
we then find, 



fii(i)/Ji = 

n a (7)M = 

nsWM = 

4 (7)M = 



M + l 

1-1/7 

1 + M7 


M|l- 7 | 
1-1/7 

f |M(2 - 7) 
\ 2 (I-I/7) 



for 7 > 
for 7 < 

for 7 > 
for 7 < 



-1/M 
-1/M 

-1/M 
-1/M 



for 1 j= 7 > -1/M 
for 7 < -1/M 

fl| for 7 > -1/M 
for 7 < -1/M, 



(31) 
(32) 

(33) 



and for antiferromagnetic couplings, we find, 

OW/17,1 _ f |1 - l/Tl for 7 > 1/M 

ni(i)/\Ji\ - | M _ 1 for7<1/M 

nn/ , n J for 7 > 1/M 

" 2(7)/|Jl1 = \ l-M7 for 7 < 1/M 

OW/17,1 _ ill-1/71 for 7 > 1/M 

V 3 (l)/\Ji\ - [ M(1 _ 7) fo r 7< i/M 

n M/\T\ - / 2 1 1 1/^1 for 7 > 1/M 

SU(7)/|Ji| - (m(2- 7 )-1 for 7 < 1/M. 



(34) 

(35) 
(36) 
(37) 

(38) 



We remark that these formulae also apply for the isosceles 
triangle, M = 2.@ 



V. LOW TEMPERATURE NUMERICAL 
RESULTS FOR N = 4 

In Fig. 1, we plot the mode frequencies ^(7) relative 
to for the squashed tetrahedron case M = 3. The 
upper and lower panels correspond to the FM and AFM 
cases, respectively. The circle in the upper panel of Fig. 
1 denotes the absence of a zero-frequency peak at all tem- 
peratures for the regular tetrahedron. We have verified 
these mode frequencies by numerical evaluation of the 
explicit integral representations of 8Cn{us) and SC^uj). 
For example, in Fig. 2 we show the low-T behavior of 
SCu(ui), presented as log 10 [5Cn(u;)] versus w/|Ji|. For 
the FM case with 7 = 0.3 at a = 50 pictured in the 
upper panel of Fig. 2, 8C\i{lo) exhibits very sharp peaks 
at the frequencies f2j, where f2j/Ji = 4,1.9,2.1, and 6.1 
for i — 1, ... ,4, respectively. ^44(0;) has a single sharp 
mode at the frequency f2j_. This figure also shows that for 
7 = 0.6, the FM 8Cn(ui) modes are also sharp at a = 50, 
appearing at 4, 2.8, 1.2, and 5.2, respectively, and at 
7 = 0.9, they appear at 4, 3.7, 0.3, and 4.3, respectively. 
We note that the mode is much weaker in intensity 
than the other modes at this temperature. For the AFM 
case, the modes tend to be much broader, as pictured for 
7 = 0.6 in the bottom panel of Fig. 2. In this case, the 
low-T mode frequencies satisfy fij/| JjJ = 2/3, 0, 2/3, 4/3, 
so that fii and r2 3 are degenerate. This degeneracy is ev- 
ident in the shape of the combined mode, which appears 
to consist of two peaks with different widths, both cen- 
tered at w/|Ji| = 2/3. In addition, is a central peak, 
which grows in intensity as T decreases. 

For the special case of the three-pronged star, 7 = 0, 
the leading behaviors of the low-T modes are presented 
in Fig. 3. For the FM star, pictured in the upper panel of 
Fig. 3, we have plotted n 2 8Cu{Lo) versus \a\(u>/Ji — n)/n 
for the largest amplitude modes Oi and f^, for i = 1,4. 
Since fix and Q2 appear at lo/ J\ = 4, 1, respectively, and 
since the VL\ modes present in 8Ca and 8C\\ are weaker 
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FIG. 1: The low-T magnon mode frequencies for the FM 
(top) and AFM (bottom) cases for the squashed tetrahedron 
(M = 3). 



FIG. 2: Plots of log 10 [(5Cii(w)] versus u/\J\\ for the squashed 
tetrahedron (M = 3) at very low T. Top: FM case at a = 50 
for 7 = 0.3, 0.6, 0.9. Bottom: AFM case for 7 = 0.6 at various 
low T values. 



than the ^2 mode, this presentation was chosen for clar- 
ity. Each of these modes was plotted at a = 5, 10, and 
20, demonstrating the low-T scaling that occurs. We also 
note that the ^2 mode in SCh(lv) drops discontinuously 
by many orders of magnitude (and to zero as T — > 0) 
at u>/Ji = 1, as indicated by the ~ sign. This behavior 
is very similar to that of the FM chain, except for the 
difference in the frequencies involved. 

The AFM three-pronged star has parameters close to 
those present in the squashed tetrahedron Fe4.0,Q The 
strongest low-T modes are pictured in the bottom panel 
of Fig . 3 , in which we plotted 8Ca (uj ) versus | a \ (ui/ \ J\ \ — 
n)/n for i = 1,4, n = 1, 2, and a = —5, —10, —20, and - 
40. Since these modes are sufficiently close in magnitude, 
the 8Cu(ui) are not scaled in this figure. The additional 
modes at u)/\ J\ \ = 3, 5 are very weak, and are not shown. 
As for the FM case, SI2 drops discontinuously by orders 



of magnitude at lo/\Ji\ = 1, vanishing as T — > 0. In ad- 
dition, in both cases, the mode shapes approach uniform 
functions of |a|(cj/| Ji| — n) as T — * 0. This behavior is 
actually simpler than that obtained for the AFM chain, 
[8| because in that case, the Qi and O2 modes present in 
SCii(ui) both approach the same frequency, w/|Ji| = 1, 
as T — > 0, making it difficult to separate them. 

It is interesting to compare these findings with the 
simpler results in the case of a perfect tetrahedron (the 
equivalent neighbor model with N = 4. There only 
one low-T mode is present, at 51/ J = 4 in the ferro- 
magnetic case, or at = in the antifcrromagnetic 
case. The low-T scaling of these single modes was shown 
previously. \u\ Allowing one spin to be coupled differently 
induces a splitting in the spectrum of low-T magnons, a 
phenomenon which was alreadyobserved in the study of 
the isosceles triangle of spins. |8( 



6 
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FIG. 3: Top: Plots for M = 3 of the fi a and ft 2 modes in 
n 2 5Ca(io) versus |a|(w/Ji — n)/n with i,n= 1,4 for the FM 
star, 7 = 0, at a = 5, 10, 20. Bottom: Plots for M = 3 
of the Qi and ^2 modes in 8Cu(u>) versus |a|(w/|Ji| — n)/n 
with i = 1,4 and n = 1,2, for the AFM star, 7 = 0, at 
a = -5,-10,-20, and -40. 

In Fig. 4, we plot the full temperature dependence of 
the two primary modes present for the AFM case with 
7 = —0.05, which is thought to be a better approxima- 
tion to the parameterspresent in Fe4 than in the bot- 
tom panel of Fig. 3. [J, |5j In this figure, we show the 
results of calculations for both SCh(lu) and SC^lo) at 
a = 0, -1, -5, and -20. At infinite T, a = 0, 5C aa (uj) 
exhibits a broad peak with a maximum at u>/\ J\\ « 1.7, 
and 5Cu(ijj) has substantial weight at low frequencies, a 
well-defined peak at lj/\J\\ « 0.8, and a small peak at 
w/\J\\ ~ 1.05. As T is lowered, the peak in SC^lo) de- 
velops into the sharp Q,\ mode, approaching oj/\ J\\ = 2 
as T — » 0. In addition, <5Cn(w) develops into the two 
modes ^2 and fli at cj/\Ji\ = 1-05 and 2, respectively. 
The minor peaks at u>/\Ji\ w 3.15 and 5.15 are too weak 
to show up on the scale used in this figure. 
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FIG. 4: Plots for M = 3 of 5Cu(ui) versus w/|Ji| at various 
temperatures, for the AFM case with 7 = —0.05, appropriate 
for Fe 4 . B 

Finally, in Fig. 5 we show low-T plots at N = 4 of 
5Cn(w) versus wlal 1 / 2 /! J±\ for the special points 7 = 
±1/3, corresponding to the onsets of the central peak of 
the mode ^2- In both cases, curves for |a| = 160, 1280 
are shown. Remarkably, the FM and AFM cases are 
nearly identical, when plotted in this manner. As for the 
similar scalings at the endpoints of the the parameter 
range of the central peak for the isosceles triangle, || 
this scaling only applies to the frequency, without a cor- 
responding scaling of SCh(lu), so that the overall scaling 
does not correspond to a scaling of the time in 8Cn(t). 
However, for the isosceles triangle, the FM and AFM 
cases appeared to be nearly similar at temperatures that 
differed by a factor of about 8, whereas for the squashed 
tetrahedron, the temperatures are essentially identical. 



VI. CONCLUSIONS 

We have solved for the time correlation functions of 
the A^-spin squashed equivalent neighbor model, with one 
spin coupled via the classical Heisenberg exchange 3\ to 
the M = N — 1 other spins, all of which are coupled to 
each other via a different Heisenberg exchange J2. Our 
results are qualitatively similar to those of the isosce- 
les triangle, N — 3, but show that for arbitrary N > 3, 
there are only four low-temperature modes, given by Eqs. 
(|31|l - H34f) and (|35|I - H38|I for ferromagnetic and antiferro- 
magnetic signs of J\ , respectively. 

At infinite T, we showed explicitly that the long- 
time asymptotic behavior of the autocorrelation function 
Cji (t) on a prong of an M-pronged star approaches 
its asymtotic limit as We also showed that for 

3 < N < 8, the infinite-T, long-time asymptotic limit 
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FIG. 5: Plots for M = 3 of 6Cu(u) versus u\a\ 1 ' 2 /\.h\ at 
\a\ — 160, 1280 for the Q2 modes at the onsets of the central 
peak, 7 = ±1/3 for the FM and AFM cases, respectively. 



of Cn n (t) is greater than that of Cji(t), and speculate 
that this relation is likely to hold for arbitrary N. We 
also showed that at infinite T, C^f 1 ^ 1 ^) approaches its 
long-time asymptotic limit as (t)~ M , one power slower 
than does Cjvjv(i)- 

We showed explicitly that these mode frequencies ap- 
ply for the isosceles triangle (N = 3) and for the squashed 
tetrahedron (N = 4). 8] For the particular parame- 
ter values appropriate for the single molecule magnet 
Fe4, with four 5 = 5/2 Fe +3 spins on the corners of a 
squashed tetrahedron, we expect that this classical cal- 
culation of the Fourier transform of the time correlation 
functions will represent a reasonably good envelope of the 
(^-functions present in the quantum mechanical treatment 
of this model, provided that the temperatures are not too 
low with respect to \Ji\. Thus, wc expect the qualitative 
features shown in Fig. 4 and the lower panel of Fig. 3 
to be observable in inelastic neutron scattering studies of 
single crystals of Fe4. 



APPENDIX 



The integrals appearing for M > 3 in Eq. (|18J) are 



h(t) 



(^lzo)i 
42 o2 



(~<1 



cos(s<*)), 



(A39) 
(A40) 



h(t) = -<(AS lz0 ) 2 cos[(l- 7 )^1}, (A41) 
cos{[s - (1 -i)x]t*} 



(x - C x ^ M f 



(A42) 



The density of states for N spins is given by 



V N {x) 



-2p(x) = ^ 



E[(N-l)/2] 

G(x) ^ <d(N-2p-x) 

p=0 

x8(i — N + 2p + 2)d N - 2p (x), 

p {-\) k {N -2k-x) N - 2 (N 

k 



(A43) 



2 N - 1 (N- 2)! 



(A44) 



where E(x) is the largest integer in x and Q(x) is the 
Heaviside step function. As noted in Eq. (|13|l . the Ii{t) 
contain integrations over T>M-i{y)- 

Although the Ii(t) are explicitly triple integrals over x, 
y, and s, the only y dependence of the integrand appears 
in the expressions for S 2 z0 and (ASi z0 ) 2 , given by Eqs. 
(|llfl and (|12|) plus the expressions following Eq. Q. In 
most of these integrals, one has to evaluate 



h = 



M-l 



dy 



v+i 



Is/ — 1 1 



dxg M (x,y)f(x,t) (A45) 



M — 2 



dxf{x, t) 



dxf(x, t) 



M-2 

where gAi(x,y) has 
1 - 2/ 2 ) 2 /(4^)]P M -i 



dygM(x,y) 

-i| 

M-l 

dyg M (x,y), 

c-l\ 



either the form a(x)[l — 
(y) or the form a(x)(x 2 



(A46) 

(x 2 + 
+ 1 - 



y ) T>M-i{y), and f(x,t) involves an integral over s. In 
most cases, the y integrals can be performed before the 
x integrals, reducing the triple integrals to double inte- 
grals, precisely as was done for the equivalent neighbor 
model with M -> N.ftt\ 

We now calculate the exact infinite-time, infinite- 
temperature limit of the correlation function Cn(t) for 
7^1 from Eq. (|A39(I . We first perform the integration 
over s, and then invert the order of the remaining two 
integrations, as outlined above. For N = 4, we then find, 



lim Clt\t) 



29 
360 

1 

^96 
1 

~192 
0.355496 



7T 

384 
Li 2 

In 



83 
360 

1 

2 
n 2 



In 2 



Li 2 



Is 1 " 3 



(A47) 



where Li2(z) is the standard dilogarithm function, 

r0 ln(l - t) 



Li 2 (z)=^ 



t 



-dt. 



(A48) 



The exact formulae become increasingly complicated 
with increasing N, so in Table I, we only list the nu- 
merical values of those additional ones for 3 < N < 11, 
along with those of the infinite t, T limits of CjvA r (0- 
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N 




lim t->oo Cjvjvfi) 

T-.00 


3 


0.370130 


0.480521 


4 


0.355496 


0.436345 


5 


0.342702 


0.416362 


6 


0.337024 


0.401888 


7 


0.333611 


0.384419 


8 


0.331595 


0.378635 


9 


0.330327 


0.374027 


10 


0.329516 


0.370270 


11 


0.328992 


0.367148 



TABLE I: Infinite t, T limits of the autocorrelation functions 



Next, we sketch our procedure for obtaining 
limT^oo SC°f 7#1 (£) as t > 1 for arbitrary N. From Eq. 
(|A41|) . we first perform the integration over s, giving us 
a function proportional to f(x) given by Eq. I|25|) . To 
avoid the singularity at x = 1, we do not invert the order 
of the remaining two integrals, but instead integrate with 
respect to x by parts twice, leading to 



lim 8Clf'\t) 



1 



E 



8t CT t±i 



M-l 



ydyV M -\{y) 



xf(y + a)cos[(y + a)t}. (A49) 



We then integrate with respect to y a total of M — 2 
times, noting that all terms proportional to derivatives 
of / sum to zero. We finally obtain, 



lim SCI? ' 1 ® 



1 



E(M/2 

E 

p=0 



A 



Mp 



(t) M 

x cos[(M 



2 M+1 

x (M - 2p 



A Mp f{M - 2p) 
-2p)t + Mn/2],(A50) 

(1 - Sp,M/2) 



V, 



M - 1 
P 



(l-<y Pl0 )(M-2p + l) 
M — 1 N 
p-1 



(A51) 



where /(x) is given by Eq. We note that /(0) = | 

and /(l) = 2. 

We now sketch our derivations of the low-temperature 
mode frequencies. We first note from Eqs. (|17fl and 
(|A40|) that the Fourier transforms of <5Cjvjv(£) and I\(t) 
both contain S(s — ui), where u> = u)/\J\\. From the above 
discussion, each of these then can be reduced to a single 
integral over x, 



K {Q) = 



min(M,tD+l) 



cLxQm{x, oj)e 



a(i-l)x 2 +ui' 2 ] 



(A52) 



where Qn(x, w) is different for 5Cnn{u) and the f2i mode 
contribution to SCh(uj). In both cases it is independent 
of a and T, and is therefore irrelevant to the determi- 
nation of the mode frequency f?i in the limit T — > 0. 
The integration limits arise from the condition that the 
(5-function is restricted by|i — l|<s<x + l. For the 
FM case, a > 0, we first consider the case 7 < 0. As 
a — > 00, the integral is maximized by choosing x to have 
its minimum value, x — \u) — 1|. We then maximize the 
resulting expression for the exponent as a function of u>, 
which occurs at ui = ui* = 1 — I/7. For 7 > 1, the 
minimum x value, \u> — 1|, is limited for large w by M, so 
Hi* = N. The crossover occurs when these frequencies are 
equal, N = I-I/7, c*7= -1/M. Settings* = fii/|Ji|, 
we thus recover Eq. (|3T|l . For the AFM case as T — > 0, 
a — > —00, we want to minimize (7 — l)x 2 + uj 2 in the 
exponent. For 7 > 1, this occurs at x = \u) — 1|, and for 
7 < 1, it occurs at x = uj + 1. In both cases, optimizing 
the exponent leads to uj* = |1 — I/7I. The latter case is 
restricted by the limitation uj* = M — 1. The crossover 
between these two limits occurs at M = |1 — I/7I, or 
7 = 1/M. Setting uj* = fii/|Ji|, we then recover Eq. 

(EJ. 

We now focus on the integral ^(t), Eq. I|A41(I . We 
first perform the y integral as sketched above. Then, the 
integral over s does not contain any time dependence, 
and as T — > 0, it is dominated by the factor exp(cvs 2 ). 
After integration by parts, we obtain the single integral 
over x, which has the form 



cIxPn(x) exp[a(7X 2 ± 2x)\ cos[(l — ^xt*], 



(A53) 

where Pn{x) is independent of a, as in Eq. (@J. Fourier 
transformation then involves the (5-function, 8(uj — |1 — 
7|a;), so that the position of the mode due to I2 is 
found by optimizing the expression exp{a[7<I> 2 /(l — 7) 2 ± 
2a)/|l — 7I]}. For the FM case and 7 < 0, we maximize 
this function with the + sign, leading to a)* = 1 — I/7. 
For 7 > 0, the ^-function was restricted by x < M, lead- 
ing to uj* = M\l - j\. These values for f2a/|Ji| = uj* 
are equal at 7 = —1/M. Combining, we obtain the FM 
f?3 mode frequencies, Eq. I|33|l . For the AFM case as 
a — > —00, we choose the — sign in the above exponent, 
and minimize 7<I> 2 /(1 — 7) 2 — 2uj in the exponent. For 
7 > 0, this occurs at uj* = |1 — For 7 < 0, the over- 
all exponent is bounded by d)*/|l — 7I < M . Combining, 
we obtain the expressions for f^/l Ji| for the AFM case, 
Eq. (EH). 

We now turn our attention to I3. In taking the Fourier 
transform, there are four (5-functions, 8[u> — s — (1 — j)x], 
5[Q+s+(l— 7)2;], 5[oj+s— (1—7)0;], and S[uj — s+(1— 7)2;]. 
These ^-functions lead after the usual reductions of the 
y integrals to the following integrals, respectively, 

|.min[Af,(<i+l)/(2-7)] 

K\{uu) = I cIxRni{x,uj) 

Jmax[(l-cD)/ 7 ,(Li-l)/(2-7)] 
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xf+(x,u>), (A54) 

,nun[Af,e(7-2)(tD+l)/(7-2)] 
K 2 (uj) = 6(7-1) / dx 
Jmax[0,(l+tD)/7,(ii-l)/(7-2)] 

xR N2 (x,Lj)f-(x,Lu), (A55) 

<-min[itf,(l— oi)/7] 

X- 3 (w) = 6(1 - 7 ) / dx 

Jmax[0,-(l+£D)/7,(£D+l)/(2-7)] 

xR N3 (x,ui)f + (x,ui), (A56) 

/■min[M,(a>+l)/7] 
-^(w) = / dx 
Jmax.[0,(Q— 

xR Ni {x,Q)f-{x,w), (A57) 
f±(x,Q) — exp{a[(i 2 + 7(7 — l)a; 2 ± 2(D(7 — l)x]}, 

(A58) 

where the RNi(x,uj) are independent of T. 

We first consider the AFM case of K\, a — > — 00. For 
7 > 1, all terms in the exponent are negative, so we 
need to minimize the function uj 2 + 7(7 — l)x 2 + 2xd>(7 — 
1). Setting x = (uj — l)/(2 — 7), and optimizing this 
function with respect to w, we find that its minimum 
occurs at uj* — 2(1 — 1/7)- For 7 < 0, the last term 
in the function to be minimized is negative, so we take 
x = (uj + l)/(2 — 7) from the upper integration limit. 
Optimizing the function, we find Q* = 2(1/7 — 1), so both 
7 regions satisfy uj* = 2 1 1 — 1/7 1 . However , this is sub j ect 
to the constraint on the upper integration cutoff, which is 
(d)* + l)/(2-7) = M, oru>* = 2M-I-M7. These values 
are equal at 7 = 1/M. Altogether, uj* — fi4/|Ji| for the 
AFM case in Eq. (jSHJ- For the FM case with 7 < 0, 
we take x — (uj + l)/(2 — 7), optimize, and again obtain 
uj* = 2(1 — I/7). The cutoff occurs when the lower limit, 
x = (£-l)/(2-7), equals M, giving uj* = 2M+I-M7. 
The crossover occurs at 7 = 1/M, as given by Eq. 
forf7 4 /|Ji|. 

Next, we consider the FM case of K4. First for a — > 00, 
7 < 0, it is easily seen that the exponent in /_(xtD) is 
positive definite. Thus, we might expect the upper limit 
for x to apply. But, this is either the cutoff, M, or a 
negative quantity, (uj + l)/7- Thus, the only positive 
limit is the lower cutoff, x = (Q — l)/7, which can be 
positive for uj < 1, leading to a larger exponent than 
obtained by setting x = 0. However, /-[(d) — 1)/j,uj] — 
exp[(a/j)(uj 2 + 7 — 1)], which for a > 0, 7 < has a 
maximum at uj* — 0, corresponding to a central peak. 
This will be the mode frequency until (uj — = M, 
the upper cutoff, resulting in uj* = I+M7. The crossover 
occurs at 7 = —1/M. Thus, this mode reduces to f2 2 /| Ji\ 
as given by Eq. (|32l . For the AFM case for 7 > 0, we 
set x = (uj + l)/7, and again we find /-[(cj + 1)/j,uj] — 
exp[— (\a\/j)(uj 2 +7 — 1)], which has a maximum at uj* — 



0. This form continues until (1— uj)/j = M, which occurs 
at uj* = I-M7. The crossover occurs at 7 = 1/M. Thus, 
this gives rise to the mode f22/|<A| in Eq. 

We now consider the Ki integral. This makes a very 
small contribution, because of the the severe limitation 
that it vanishes unless 7 > 1. For the AFM case, the 
exponent is optimized at x — x* = tD/7, and then op- 
timizing the mode frequency with respect to a), we find 
that uj* = 0, so that K2 for AFM coupling contributes 
to f^/l Ji|. For the FM case, the maximum exponent oc- 
curs at x = M, and from the (5-function restrictions, we 
see that K2 makes a contribution to the O4 / 1 «/i | mode. 

Finally, we discuss briefly the K% case, for which 7 < 1. 
Setting 7 < for the FM case, the optimum situa- 
tion is obtained when uj* — 0, so that it adds to the 
VI2/ 1 Ji I mode. For the AFM case, the optimum x value 
is x* = —uj/j, and this is restricted by x < M. Hence, 
K3 essentially makes a contribution to the Q2 mode for 
the AMF well. 
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